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We study the high latitude {\b\ > 10°) diffuse 7-ray emission in the Galaxy in light of the recently 
published data from the Fermi collaboration at energies between ICQ MeV and 100 GeV. The 
unprecedented accuracy in these measurements allows to probe and constrain the properties of 
sources and propagation of cosmic rays (CRs) in the Galaxy, as well as confirming conventional 
assumptions made on the interstellar medium (ISM). Using the publicly available DRAGON code, 
that has been shown to reproduce local measurements of CRs, we study assumptions made in 
the literature on atomic (HI) and molecular hydrogen (112) gas distributions in the ISM, and non 
spatially uniform models of diffusion in the Galaxy. By performing a combined analysis of CR and 
7-ray spectra, we derive constraints on the properties of the ISM gas distribution and the vertical 
scale height of galactic CR diffusion, which may have implications also on indirect Dark Matter 
detection. We also discuss some of the possible interpretations of the break at high rigidity in CR 
protons and helium spectra, recently observed by PAMELA and their impact on 7-rays. 
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I. INTRODUCTION 

The study of the physics of Galactic cosmic rays (CRs) is one of the most active research areas at present. Sensible 
advances in the field have come in connection to the wealth of high-accuracy data recently collected by several new 
instruments, with further progresses expected in the upcoming future. In particular, since its launch three years ago, 
the Fermi Gamma-Ray Telescope [l| has been producing the most detailed and precise maps of the 7-ray sky ever, 
given its wide energy coverage and excellent energy resolution, its large effective area and field of view, as well as the 
best angular resolution for a 7-ray detector in space (for details on the performances of the instrument see 0). Since 
the interaction of Galactic CRs with the interstellar medium give rise to a copious 7-ray yield, the diffuse emission 
in the Milky Way is by far the brightest source detected by Fermi. It is then expected that Fermi data will drive a 
significant improvement in understanding the origin and propagation of CRs. 

A key feature that would be particularly important to establish is whether there is room for (or eventually even 
need for) an exotic source of 7-rays and/or CR leptons in the Milky Way, on top of the astrophysical sources most 
plausibly providing the bulk of Galactic CRs, namely supernova remnants (SNRs) and, possibly, pulsars. Such exotic 
component is predicted, e.g., in connection to the dark matter halo of the Galaxy for several dark matter candidates, 
the prime example being Weakly Interacting Massive Particles (WIMPs), i.e. early Universe relics which have a small 
but finite probability of annihilating in pairs and produce multi-GeV (or -TeV) yields (for a recent review on WIMP 
dark matter, see, e.g. @; the issue of probing and constraining dark matter models with Fermi-LAT data was recently 
discussed, e.g., |4l-lll|). Suggestions of an improvement in the fit of Fermi 7-ray data at intermediate latitudes taking 
into account the spectral hardening due to a signal from annihilating dark matter have already been proposed in the 
literature, see, e.g., [l^l- On the other hand it is evident that for such kind of analysis it is crucial a very accurate 
modeling of background components, whether of galactic, extragalactic, or instrumental origin. 

The question we wish to address in this work is to what extent the cross correlation of the 7-ray data with other 
local CR probes, namely the local measurements of the flux of CR primaries, the ratio of secondary to primary CR 
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nuclei and the flux of electrons and positrons, helps disentangle degeneracies between the type and distribution of 
sources and the propagation modeling similarly to what was done in the past with the EGRET data by The 
emphasis on the issue of locality comes from the observation that, in the vast majority of works in the literature, 
propagation of CRs in the Galaxy is treated with an effective approach in terms of a diffusive/convective equation 
encoding, on average, through a set of simplifying assumptions and few parameters to be tuned to the data, the 
physical process of charged particles making a random walk in the regular and turbulent Galactic magnetic fields. 
Local measurements give fairly good tests of average properties of the local medium, with the physical averaging scale 
depending on the species considered, and with the caveat that some of the parameters in the propagation model show 
patterns of degeneracy, as studied at length in the literature. After selecting models fulfilling these local constraints, 
covering a wide range of different physical propagation regimes, as well as gas distributions and, to some extent, 
galactic (but excluding dark matter) source distributions, we will derive predictions for the 7-ray flux at intermediate 
and high latitudes, compare against Fermi data and discuss whether a discrimination among the different models 
is possible. The choice of excluding from the analysis low latitude data is connected again to the issue of focussing 
on local properties, since although 7-ray data naturally reflect a global observable, summing along the line of sight 
contributions to the emissivity from all regions of the Galaxy, high latitude fluxes are dominated by the local terms. 

The Fermi Collaboration has first published spectral data at intermediate latitudes in [Tsj disproving the presence 
of a GeV excess in the diffuse 7-ray spectra suggested by EGRET [l^ . Wc will consider in the following the updated 
published spectral data from the Fermi Collaboration at latitudes 10° <| 6 |< 20°, 20° <| 6 |< 60°, | 6 |> 60° and 
0° < / < 360°, and energies up to ^ 100 GeV Such spectral data have been evaluated using a more clean 

sample of 7-rays than the publicly available "Data Clean 4" (P6V11) [l^; furthermore the CR contamination to that 
remaining data-set has been modeled in and its contribution (isotropic) to the spectra has been subtracted, thus 
making these spectra ideal for studying the large scale properties of the Galaxy. For E-y > 100 GeV, statistical errors 
and contamination of CR electrons and nuclei result in great uncertainty on the exact spectrum of the 7-rays p7| (see 
also [H]). 

Since we do not include the data from \ b \< 10° in our analysis, wc mainly probe the properties of CR propagation 
and to some extent the gasscs but we are less sensitive to the properties of the distributions of galactic sources. The 
assumptions made on how CRs diffuse away from the disk, can have a significant effect on the 7-ray distribution 
from inverse Compton scattering (ICS) by electrons/positrons that are either diffusive shock accelerated ISM at 
SNRs, secondaries from inelastic pp and pHe collisions (predominantly), or produced in Pulsars magnetospheres and 
accelerated within Pulsar Wind Nebulae (PWNe) . Also since magnetic fields decrease as we move radially away from 
the galactic center, we expect a radial profile as well in the diffusion |23| of CRs which can also have an effect on 
the observed 7-ray spectra. Gamma-rays from tt*^ decays produced in inelastic collisions of CR protons and heavier 
nuclei with ISM nuclei, and bremsstrahlung off electrons are dominant contributions to the 7-ray spectrum from 100 
MeV up to at least 50 GeV. Thus the observed 7-ray spectra can also be used to probe the ISM properties, such as 
confirming models for the HI and 112 gas distributions. 

Studies of the sky observed by Fermi, at latitudes |^> 5° have shown an excess of diffuse 7-rays towards the 
center of the Galaxy and up to latitudes of | 6 |« 50° l2l] . [22| known as the "Fermi haze" or "Fermi bubbles" that 
can be of either astrophysical [22l - [25| or DM origin [l^, l26l |. Even though this feature(s) extends up to high latitudes, 
since it is confined in longitude within | / |< 20° its effect on the spectra in the regions of our study is washed out. 
Thus our analysis can not probe the properties of the Galaxy in that region. Nor it can in the region of (| |< 10°), 
where [27l - l29j have suggested the presence of a DM signal. 

This paper is organized as follows. In section [III we summarize the assumptions that we make on the primary CR 
sources, diffusion, magnetic fields and gas models, as well as briefly present the tool we use to solve numerically the 
CR propagation equation, namely the DRAGON code [s^l . In section IIIII we describe the analysis that we conduct 
in constraining the CR propagation and ISM properties, by fitting to the local fiuxes of CRs. Section ITVl is devoted 
to studying the effects that the various assumptions on diffusion of CRs, ISM gasses, distribution of CR sources and 
7-ray production cross-section from pp-collisions have on the 7-ray spectra. We also place constraints imposed on CR 
propagation models and on models of ISM gasscs distributions using the diffuse 7-ray data[l3|, as well as discuss the 
implications of our findings. Finally in section |V] we present our conclusions. 
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II. COSMIC RAY PROPAGATION 



As suggested in [3l[ the propagation of CRs in the Galaxy at energies below 10^'' eV can be described by: 
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where il){r^p,t) is the CR density per unit particle momentum, or in terms of phase space density /(p), 'ip{p)dp = 
Anp^ f{p)dp, q{r,p,t) is the source term including primary, spallation and decay of heavier CR species. Dxx{^ is 
the diffusion tensor in physical space and Dpp{f) the diffusion coefhcient in momentum space, p is the momentum 
loss rate due to interactions with ISM, the Galactic magnetic field or the interstellar radiation field (ISRF), V is the 
convection velocity, and rjrag and T^ecay are the time scales for fragmentation loss and radioactive decay respectively. 

For our simulations we use DRAGON [sO, Is^ - ls^ ] that numerically solves Eq. [T] in the steady state approximation 
dip /dt = 0, assuming cylindrical symmetry, in a 2+1-D grid where each point is described by its galactocentric radial 
distance r g (0, 20) kpc, the distance from the galactic plane z S (— i, L) with L < 20 kpc and momentum p. 



A. Primary Sources of CRs 

We consider SNRs as the primary CR sources up to energies of ~100TeV [l^. For each nucleus i of charge Z, the 
source term describing the injection of CRs in the ISM is given as a function of rigidity, R = p/Z by: 



qi{r,z,E) = fsir,z)qo, 



(2) 



where qo.i is the normalization of the injected CR species, and fs{r,z) traces the spatial distribution of SNRs. In 
our reference model we consider the distribution by [36| , which is derived on the basis of pulsar and progenitor star 
surveys. Alternative models will be considered as well (see section IIVE|) . 

Being interested in high energy 7-ray fluxes, we pay particular attention to the proton spectrum, given that protons 
provide the dominant contribution to the diffuse 7-ray spectra in the whole energy range we consider fsee III C]) . We 
allow for a CR proton injection spectrum described by a broken power-law: 
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with two breaks at i? = i?o 1 ~ 10 — 30 GV and R ~ R^ 2 ^ 300 GV with the spectral index 7^ lying in the range 
1.85 — 2.1 at low rigidities, 7^ in the range 2.3 — 2.5 at intermediate rigidities and 73 ranging in 2.18 — 2.35 at high 
rigidities. This choice is motivated a posteriori by fitting the local proton data and in par ticular from the combined 
fit of the recent PAMELA proton spectrum [s^l, and CREAM spectrum above 2.5 TV j38|. 

Electrons and positrons accelerated between a _2ulsar and the termination shock oi^ the wind nebula, may also 
contribute to the high energy spectrum [3/ 



pulsars were found to be particularly well suited 




and then to the 7-ray flux [43|,|4J]. In particular, middle aged 
45j . Each pulsar contribution to the fluxes can be described 
by an injection spectrum ~ E~" with a high energy break Eh which is estimated at the time the surrounding PWN 
is disrupted leading to the escaping into the ISM ^ Furthermore, each pulsar has an initial rotational energy 
Wo of which only a portion 77 is injected into the ISM as CR e^. The ranges for those parameters within different 
pulsars are very broad. Indeed, n can range between 1 and 2, Wq ~ 10^^~^° erg, 77 ~ 0.1 [H, and Et ^ 10 TeV. 
However, the actual observed flux of in our position from a pulsar has a break that is related to the cooling time 
(from ICS and synchrotron radiation) of the during their propagation in the ISM [s^, 14^. To account for these 
effects, we choose to fit the properties of a pulsar distribution following the parametrization of [42j : 



Qp{r,z,t,E) = JoE--e-''/^'fpir,z). 



(4) 



^ Even though it is clear that the higher energy escape earher into the ISM, the differences in the estimated time scales for of 
energies between 1-10* GeV to escape into the ISM are negligible compared to the propagation time from the PWN to us [43. l46l . 
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M is a "statistical" cut-off, n the injection index for the distribution of pulsars,^ 



Jo 



vWoNb 



r(2 - n)A'p-^-Vgai 
(see Eq. 24 of [42|) with TVf, the pulsar birth rate in the Galaxy and 
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fp{r, z) describes the spatial distribution of young and middle aged pulsars. Since pulsars have typical kick 3D speeds 
of 200-400 km/s [13,143 , a 10^ yr old (middle aged) pulsar would move away from its original position by ~ 30 pc, 
and thus the spatial distribution of middle aged pulsars is practically identical to that of their birth distribution in 
the Galaxy as given in (ioj : 
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with Ri = 0.55 kpc, zi = 0.1 kpc, a = 1.64 and 6 = 4. 
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B. Diffusion and Magnetic Fields 

Our galaxy is permeated by a large scale, so called regular, magnetic field, and by a randomly varying, so called 
turbulent, magnetic field with comparable strength on the disk. The large scale galactic magnetic field is generally 
assumed to be a bi-symmetrical spiral with a small pitch angle [50l [. Here we assume that the regular magnetic field 
is purely azimuthal, Bq = Bqc/), and has the form 

Bo = Bh exp exp (^-^^ . (8) 

Based on the analysis of WMAP synchrotron intensity and polarization data in [5l|, as well as works including 
extragalactic rotation measures fsol. Ts^ - ls^ l . we choose Bh = 3 /xG and r^i = 11 kpc, with vertical scale Zh = 2 kpc. 
Although these values are affected by large uncertainties, they have little impact on our analysis, since the magnetic 
field enters only in the electron energy losses, which are anyway dominated by the ICS losses above few GeV, as we 
show in Appendix [XI 

The diffusion tensor can be in general decomposed in a component parallel to the direction of the regular magnetic 
field, and a component describing diffusion perpendicular to the regular magnetic field, D±. It can be shown (ssj 
that, assuming cylindrical symmetry and that the regular magnetic field is azimuthally symmetric, parallel diffusion 
is irrelevant and only D± has an effect. We consider then for simplicity that the diffusion is described by just one 
quantity, the diffusion coefficient. 

The diffusion coefficient is in general expected to depend on the position, because turbulence is not uniformly 
distributed in the Galaxy. 

In a phenomenological approach, we choose D to be described by: 

D{r, z, R) = I^o/?" (1^) ' exp (^^^) exp (^^) (9) 

with the radial and vertical scales, and z^, defining the diffusion profile in the Galaxy. Such a parametrization of the 
diffusion coefficient can be motivated on large scales since the Galactic magnetic field decreases away from the galactic 



^ Due to energy losses, the fluxes of CR with E ~ TeV, at any position of the Galaxy, are typically dominated by the contribution 
from the sources within ~ O(IO^) pc from each position. Thus the CR spectra at ~ TeV energies will differ significantly between 
different positions of the Galaxy (even at the same galacto-ccntric distances). Since the IC 7-ray spectra are studied in wide regions of 
the sky (and thus the Galaxy), we care for the averaged flux, as measured by many different observers in the Galaxy. For such a 
(statistically) averaged flux from pulsars, we follow [43l where the statistical cut-off M and injection index n refer to the statistically 
averaged values. 
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center (both in r and z). Particles gyrating along ordered field lines, may scatter due to magnetic irregularities. As 
the magnetic fields become weaker, assuming that the particles gyroradius remains small enough that the particles 
probe the ordered field component, the diffusion length (and coefficient) will increase. We choose Rq ~ 3 GV as the 
reference rigidity, while S is the diffusion spectral index which is related to the ISM turbulence power-spectrum. The 
dependence of diffusion on the particle velocity, /3 = Vp/c, is naturally expected to be linear (77 = 1), however the 
analysis by [scj shows an increase in diffusion at low energies. To account for such a possibility, the parameter 77 has 
been introduced (see e.g. [s^ls^jlll])- We will also consider the case where there is a break in the diffusion spectral 
index 6. 

In addition to spatial diffusion, the scattering of CRs on randomly moving magneto-hydro-dynamical (MHD) waves 
leads to diffusion in momentum space which results in stochastic acceleration of CRs. The corresponding diffusion 
coefficient in momentum space is related to the diffusion coefficient in physical space by Dpp oc p^va/ Dxx, where va 
is the Alfven velocity, associated to the propagation of MHD waves 

C. Energy Losses and Diffuse Gamma- Rays 

CR nuclei and CR electrons and positrons lose energy during propagation in the ISM. Depending on the energy, 
electrons and positrons energy losses are dominated by inverse Compton scattering of low energy photons of the 
interstellar radiation field (ISRF) , for which we use the model of [s^ , and, to a less extent, by synchrotron radiation. At 
energies < 1 GeV bremsstrahlung, ionization and Coulomb losses become relevant. For protons and nuclei ionization 
and Coulomb losses are the dominant continuous energy loss mechanisms. 

In addition to ionization and Coulomb losses, collisions of heavy nuclei with hydrogen or helium of the ISM gas can 
lead to inelastic scattering that can also cause the fragmentation of the parent CRs. 

CRs can also undergo radioactive decays, with the radioactive isotopes created both by fragmentation of heavier 
nuclei (e.g. -'^"Be is created from B, C, N, O) and directly in CR sources (e.g. ^^Al). All these processes are incorporated 
in DRAGON 3. 

For > 10 MeV, there are three processes that contribute mainly to the diffuse galactic component. Inelastic 
pp collisions producing tt'^s which subsequently decay to 2 photons, constitute the main contribution to the diffuse 
gamma ray fiux from the Milky Way in the intermediate GeV range, and trace the ISM target distribution dominated 
by the HI and H2 gasses. The spectral shape of these gamma rays is essentially determined by the spectral shape of 
the proton spectrum along the line of sight. CR electrons may also produce 7-rays via bremsstrahlung in the ISM 
gas, or by up-scattering low energy photons [Tsl. l6TI - [65j . Since at distances far from the galactic disk the optical and 
the IR (mainly emitted by dust) photon densities are less than those close to the disk, IC 7-rays at high latitudes 
are mainly due to up-scattered CMB, and thus due to their isotropy, they set a good probe to study the CR e~ far 
from the disk. Apart from galactic diffuse 7-rays the observed fluxes include the extragalactic flux modeled in 
and galactic point sources The CRs misidentified as 7-ray events have instead already been subtracted in the 
spectral data that we use, with the remaining CR contamination being negligible up to E-y = 100 GeV. 

D. Interstellar Gas 

The interstellar matter is made up of gas and dust with an average mass ratio of 100:1 [l^. Interstellar gas is 
composed of hydrogen, helium and small contributions from heavier elements, with hydrogen observed in atomic (HI), 
molecular (H2) and ionized (HII) states. 

The three dimensional distribution of HI gas can be derived from Lyman-a, from 21-cm spectra information and 
from rotation curves [67l - [69l |. with the 21-cm line emission being due to the transition between the atomic hydrogen 
ground state levels split by the hyperfine structure. 

We will use as a reference HI gas model the recent result obtained by [zOl , but also refer to [7l|, [t^ which has been 
widely used in the literature. 

Molecular hydrogen can exist only in dark cool clouds where it is protected against the ionizing stellar ultraviolet 
radiation. It can be traced with the A = 2.6 mm (J = 1 ^> 0) emission line of CO, since collisions between the CO 



^ Given similar assumptions, for the propagation of CRs in the Galaxy, we have checked that the output of DRAGON and Galprop 16011 
agree within an accuracy of 5 — 10% in the total galactic diffuse 7-ray flux; and by an accuracy of no more than 5% in the total 7-ray 
flux, at energies between 100 MeV and 300 GeV in the three windows that we study (0° < Z < 360°, 10° <| 6 |< 20° / 20° <| 6 |< 60° 
/ \ b\> 60°). 
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and H2 molecules in the clouds are responsible for the excitation of CO. The CO to H2 conversion factor, Xco which 
relates the H2 column density, Nh2 , to the velocity-integrated intensity of the CO line, has considerable uncertainties. 

We use as our reference H2 model the map provided by [t^ , assuming the conversion factor to vary exponentially 
with galacto-centric radius as: 



Xco{r) = 1.4cxp(r/llkpc) x 10^°i?2cm 



'^K^^km^^s, 



(10) 



comparing it with the models developed by |74| . 

The ionized hydrogen is concentrated in the vicinity of young O and B stars, with the ultraviolet radiation from 
these stars ionizing the ISM. It is known that the contribution of the HII gas to the total mass in the ISM is negligible 
[rsf . while its distribution is very similar to that of the free electrons in the galaxy. Thus we choose not to vary 
the averaged large scale distribution of this gas component, for which we use the parameterization of [t^, that was 
calculated for the Galactic distribution of free electrons: 



nniiir, z) = ne(r, z) = {11^)1 exp 



l£l 



(ne}2exp 
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Z2 



r2 



(11) 



We used the mean values for zi = 1 kpc, Z2 ~ 0.15 kpc, r2 — 2 kpc, r' ~ A kpc, {71^)1 = 0.025 cm~^, {11^)2 — 0.2 cm~^ 
and the minimum value for ri = 20 kpc. 

In general, the use of 2D, spatially smoothed gas distributions would not be accurate enough to interpret 7-ray sky 
maps with the angular resolution of the Fermi instrument. However, we will compare 7-ray spectra with observed 
spectra in very wide, longitudinally averaged, intermediate and high latitude portions of the sky, where the small 
scale features of the gas present in detailed 3D models are washed out on average mainly within equal latitude 
regions. Indeed, using our reference model, we checked that passing from a 2D to a 3D gas model changes our 7r° and 
bremsstrahlung results by no more than 10% in the regions of interest, as we show in Appendix [BJ with the difference 
on the predicted total galactic diffuse model being at the 5% level. Given that in deriving our physical conclusions 
we use the analysis carried out between the total 7-ray Fermi fluxes and the total 7-ray predicted fluxes which 
include the extra galactic background (EGB) and point sources, the impact of using a 2D ISM gas model is minimal. 

Finally, Helium appears to follow the hydrogen distribution with a factor He/H = 0.10 ± 0.08. Following [t^] we 
adopt a value of He/H = 0.11, which is widely used in the literature and neglect heavier nuclear species. 



III. METHODOLOGY 



As just illustrated, there are many unknowns involved in the modeling of both the ISM and the propagation of CRs. 
We are then forced to focus our discussion introducing a few benchmark scenarios for both aspects of the problem. 

One important parameter is the spectral index of diffusion, (5, inferred from the spectral slope of the secondary 
to primary ratios at high energy. It ranges between about 6 — 0.3 up to about 0.7. We will discuss mainly the 
theoretically motivated frameworks of the Kraichnan turbulence spectrum [tsI [79j corresponding to S = 0.5, and the 
Kolmogorov [sO] corresponding to (5 = 0.33. 

We choose a range for the vertical and radial scales of the diffusion coefficient, Zd and in Eq. (5] The discrete 
values that we choose in our analysis are Zd = (1,4,10) kpc and = (5,10,20) kpc, allowing for cases where 
the diffusion is highly homogeneous within our simulation volume (D ^ e'^l/^" kpcg(r-r0)/2o kpc-j ^^-^^ cases such as 
D g|z|/ikpcg(r-ro)/5kpc .^^jj^gj-g diffusiou properties vary significantly within the Galaxy. 

For most of the discussion we will neglect convective effects. However, we will also consider the effects of strong 
convective winds introducing in one benchmark model a convective velocity, directed only along the vertical direction 
outwards from the galactic plane Vc{z) = dvc/dz ■ \z\, with dvc/dz = 50 km/s/kpc. 

In the first nine lines of Tab. U we summarize the models we consider for the diffusion properties of the Galaxy. In 
the other models we instead use our reference model for propagation, but we vary the properties of the ISM gas and 
SNRs. 

Models are labeled in the following way: models with "KRAzd-rd" correspond to Kraichnan-like turbulence (S ~ 0.5) 
and fixed values of Zd and r^. In the same way, the "KOL4-20" model corresponds to 5 = 0.33, Zd ~ 4: kpc and 
rd = 20 kpc and the "CON4-20" has a significant convective velocity. The models labeled by "NS low" and "NS high" 
have the same (J, z^, rd) as our reference model, but use different H2 gas distributions than our reference to probe 
the uncertainties derived from [t^. The "Bronf model assumes an HI (H2) gas distributions modeled by (tiI. [72|. It^ . 
The "Source B" and "Source C" scenarios are used to study the effects of different source distributions in the Galactic 
disk. 

For each model with a set of values of S, Zd, rd and dvc / dz we derive the other propagation parameters {Dq, rj, va) 
by minimizing the oi B/C data, thus fitting our galactic (global) models to the local data, as we show in Fig. [1] 
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27.2 


0.755 


1.85/2.48/2.19 


12.3 


1.6/2.62 


0.44 


0.61 


0.44 


0.82 


1.05 


2.18/0.98/0.25 


1.09 


Scenario B 


0.5/0.33 


4 


20 





2.49 


19.5 


-0.363 


2.06/2.35/2.35 


14.9 


1.6/2.62 


0.67 


0.34 


0.6 


0.4 


0.74 


1.04/0.43/0.59 


0.71 


NS low 


0.5 


4 


20 





1.94 


14.6 


-0.324 


2.07/2.35/2.18 


15.4 


1.6/2.62 


0.68 


0.32 


0.44 


0.42 


0.67 


3.69/1.77/0.21 


1.34 


NS high 


0.5 


4 


20 





3.04 


24.4 


-0.411 


2.06/2.35/2.18 


17 


1.6/2.62 


0.74 


0.31 


0.57 


0.55 


0.85 


0.52/0.55/1.96 


0.94 


Bronf 


0.5 


4 


20 





3.39 


26.5 


-0.526 


2.08/2.35/2.18 


17.6 


1.6/2.62 


0.79 


0.38 


0.58 


0.52 


0.69 


1.45/1.42/3.44 


1.47 


Source B 


0.5 


4 


20 





2.49 


19.6 


-0.355 


2.04/2.34/2.18 


15.9 


1.6/2.62 


0.77 


0.3 


0.36 


0.39 


0.86 


1.26/0.50/0.43 


0.79 


Source C 


0.5 


4 


20 





2.33 


19.2 


-0.44 


2.05/2.34/2.18 


16.6 


1.6/2.62 


0.84 


0.33 


0.58 


0.53 


0.73 


0.73/0.34/1.44 


0.79 



TABLE I: The parameters for the benchmark models for propagation used for 7-ray predictions in Fig. [2j |4j O [6l [T] [12] and 1141 7'^'s are the injection indices for 
primary electrons below and above a break at 5 GV. Our reference model "KRA4-20" is also referred to in the text as "Scenario A", "NS mean" and "Source A". For 
"CON4-20" we assumed dvc/dz = 50 kms~^ kpc~^. The "Scenario B" model accounts for a possible break in the spectral index. In "NS low" and "NS high" we use 
H2 gas densities, respectively, 1 cr lower and 1 a higher than mean values of [73 |. however they share the same HI gas distribution of [tO]. In "Bronf we use the model 
of [t^I for H2 and the model of [tH . It^ I for HI gas distributions. The "Source B" and "Source C" models refer to different assumptions for the primary CR source 
distributions. See text for more details on the definition of the other parameters, x^'s refer either to the goodness of our fits of CR nuclei, protons and leptons or show 
the level of agreement of our predictions with 7-ray and antiproton data. 





Benchmark 


Fitted 


Fitted 


Fitted 


Predicted 


Name 


S 


Zd 

kpc 


Td 

kpc 


Do X 10^® 
cm^/s 


VA 

km/s 


V 


7f/7f/7f 


^0,1 

GV 


7i77i 


xlO"*^ erg 


Xb/c 


xl 


2 

X(e-+e+) 


xl 






KRAl-20 NS low 


0.5 


1 


20 


0.414 


11.9 


-0.454 


2.07/2.35/2.18 


16.6 


1.6/2.62 


0.31 


0.44 


0.56 


0.67 


0.74 


6.77/3.89/0.77 


2.43 


KRAlO-20 NS low 


0.5 


10 


20 


3.35 


14.4 


-0.331 


2.06/2.35/2.18 


15.2 


1.6/2.62 


0.95 


0.36 


0.85 


0.31 


0.68 


2.84/1.28/0.14 


1.09 


KRAl-20 Bronf 


0.5 


1 


20 


0.831 


25.4 


-0.563 


2.07/2.35/2.18 


18.3 


1.6/2.62 


0.33 


0.29 


0.9 


0.47 


0.9 


0.66/0.47/1.11 


0.81 


KRAlO-20 Bronf 


0.5 


10 


20 


5.6 


24.6 


-0.576 


2.09/2.36/2.18 


17.6 


1.6/2.62 


1.15 


0.35 


0.75 


0.58 


0.68 


1.43/1.29/3.23 


1.4 



TABLE II: Parameters for models (not shown in any Figure), that represent two extreme cases of either a thin diffusion halo with a low ISM gas density assumption: 
"KRAl-20 NS low", a thick diffusion halo with a high ISM gas density assumption: "KRAlO-20 Bronf; and two intermediate cases: "KRAlO-20 NS low" (thick 
diffusion halo with low ISM gas), "KRAl-20 Bronf (thin diffusion halo with high ISM gas). As can be seen from the "predicted" columns the 7-ray spectra are 
more sensitive than the antiproton spectra in discriminating among some of these cases. 
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FIG. 1: Reference Astrophysical modeL We assume that the difTusion coefficient D = 2.49 x 

10^** cm^s /3^°-^^^(i?/3 GV)°-^e''=l'"' kpcg(r-rQ)/20 kpc (^here R is the rigidity) as in Eq. H Values of of the refer- 
ence model for each observable are also shown. Upper left: fit to the B/C data. Upper right: proton spectrum, where we fit 
the injection spectrum using PAMELA [sj and CREAM data [3^ (see section FlI A|) . Lower left: the predicted antiproton 
spectrum, which provides a good fit to the PAMELA data [stI ]. Lower left: e~ + e"'" fiux, including Fermi observations [s^ . 
Constraints on the primary, secondary and pulsar fluxes properties are obtained by fitting to the data. Dotted lines refer to 
unmodulated CR fluxes. 



(upper left panel). We use the HEAO-3 [8l|, CRN [H and CREAM ^ data points. We then fix the spectral indices 
for protons 7^, 7^,7^ below and above the rigidity breaks Rq ^,Rq2 by fitting to the recently released PAMELA [s^l 
and CREAM data [3^ (see Fig. [T] upper right panel) . We then refit the normalization for the proton spectral data 
and modulation potential (in the force field approximation [sj]) checking also for consistency with the BESS data of 
years 1997-1999 [8^, where the only free parameter in fitting the entire spectra of different years is the modulation 
potential. We fit in the same way the injection parameters of He nuclei (not shown in Tab. |T| by fitting to the most 
recent data [13, up to the highest energies. Having reproduced primary protons and He, we check also whether 
the predicted antiproton flux is consistent with local data (Fig. [1] lower left panel). The corresponding x^'s are listed 
in Tab. HI 

Having chosen the ISM gas models and fixed the diffusion and re-acceleration properties in the ISM, as well as the 
CR nuclei spectra, the remaining task before calculating 7-ray spectra is to fix CR electrons (and positrons) source 
properties. Since the e+ + e~ spectrum below i? ~ 30 GcV is dominated by shock accelerated electrons in SNRs 
(primaries), and by secondary electrons (and positrons) from inelastic collisions of CR nuclei with the ISM, we fit 
the primary and secondary electron spectral properties to the low energy + e~ spectrum between 7-30 GeV as 
measured by Fermi [11]. Pulsars within ^ 3 kpc can contribute to the e"*" -|-e~ spectrum up to 0(0. 1) at ii^ w 50 GeV 
and up to 0(1) at i? w 500 GeV ([i^l and also similar works [10, IHl, IH, [s^ ) . Thus, if we assume pulsars contribute 
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maximally we find from the Fermi data the injection index n for the distribution of pulsars of Eq. |3]and the averaged 
total energy injected into the ISM through CR per pulsar rjWo- Best fit values are found to be n ~ 1.4 and 
M ^1.2 TeV. Using them, and assuming a constant birth rate of iVf, = 30 yr~^, wc find that on average wWn — 10^^ 
erg for our various propagation models (see Table U) , which is well within the allowed range of values ^2|, l90l - [9^ . 

Having fixed all the properties of the CR electrons from SNRs, pulsars and inelastic collisions, we can then compute 
the 7-ray diffuse spectra via up-scattering of ISRF and CMB photons by the CR , bremsstrahlung in the ISM gas 
by both CR electrons and protons and by decays of 7r° produced in pp collisions in the ISM gas. We note that while 
it is clear that the high energy part of the e~ + e+ flux is dominated by local sources, in our approach wc fit the 
properties of the statistically averaged flux [n, M and 77 Wq) to the + e+ data, and extend those properties to 
the entire distribution of [43 of pulsars in the Galaxy. If for some reason the local e~ + flux above 100 GeV is 
severely enhanced (suppressed) by the presence (absence) of strong local sources, this should have an effect in our 
model over-predicting (under-predicting) the ICS components of the spectra mainly at lower latitudes. 

We also note that we have checked for consistency with the PAMELA positron fraction [ol, [oj] and the recently 
released electron spectrum [95| . 

We bring attention to the fact that for each model we use exactly the same gas distribution model to first fit 
propagation and injection properties against primary CRs, then to predict secondary antiprotons and Icptons and 
finally to produce 7-ray maps. 



IV. RESULTS 



A. Reference Model 



As it is clear from Fig. [2l our reference model "KRA4-20" (see Table |I| provides a very good combined fit of the 
local CRs (see Fig.[T|) and the 7-rays at intermediate and high latitudes. In Fig. [3] we also show the SNRs and pulsars 
contributions to the total spectra separately. 

The best fit to the 7-ray spectra is achieved at the higher latitudes (| & |> 20°) that are also less affected by 
uncertainties in the sources distributions. This can also be seen by comparing in Fig. [3] our predictions between 
10° < I 6 I < 20° and | & | > 60° where we show separately the contribution from SNRs and pulsars that have different 
source functions, given in eq. [5] and 01 At intermediate latitudes it slightly under-predicts the observed flux while still 
giving a fairly good fit. We note that the "source" component that we show is composed by the sources detected with 
at least 14a- and also weaker sources that have been catalogued by LAT [l3|. Yet, very dim 7-ray sources that would 
be contributing, per energy bin and pixel, less photons than the uncertainty of the true diffuse background arc not 
included in the "sources" component. Such a class of sources could be millisecond pulsars (MSPs) in the Galactic 
Ridge and halo that are not accounted for. 

MSPs that are not in globular clusters, can contribute in the lower latitudes and could possibly compensate for 
our under-prediction of the total gamma-ray flux at ^ few GeV. The uncertainties in the contribution of dim MSPs 
to the diffuse galactic flux have been shown to be very significant [9^ due to the great energy loss time scale (■^lO 
Gyr [43, [93, Isl]) of MSPs which results in thei r tot al population being greatly affected by the uncertainties in the 
evolution of the Galactic halo [1^ . Recently, [9^ Il00l | have suggested that MSPs could be contributing to the isotropic 
diffuse 7-ray flux. Since MSPs occur in regions of high stellar densities such as the Galactic Ridge, and possibly the 
Galactic halo (9^ at its earlier stages, it is unlikely that the main part of their diffuse contribution is g oing to be 
isotropic. While individual MSP spectra may vary significantly, based on the measured spectra of 8 MSPs |10ll | their 
distribution spectrum could be described by 

^ ^ ^;-I^e-^/^^ (12) 

with r = 1.5 ± 0.4 and Ec = 2.8 ± 1.9 GeV and luminosity in 7-rays of L = lO^^-^io.e ei.g/s_ gy comparing to the 
data, we find that to fit to the 10° <| 6 |< 20° data we need a flux of 

w 10"^ McVcm"2s"^sr"\ between 1 and 10 GeV, (13) 

dE 

which would lead to a population of ~ 10^ MSPs in that region following the assumptions of [96| . 
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FIG. 2: Gamma-ray spectra for given sky regions predicted in our Reference Astrophysical model "KRA4-20". Upper left: 
10° <| 6 |< 20° and Q° <l< 360°, Upper right: 20° <| 6 |< 60° and Q° < I < 360°, Lower: 60° <| & |< 90° and Q° <l< 360°. 



B. Varying the Diffusion of CRs in the ISM 

In Fig. |4]we compare the 7-ray spectra predicted by models with different diffusion spectral index 5 in the three 
sky regions under study. Changing 5 affects the proton spectra whose propagation timcscalc depends on the diffusion 
timcscalc. Lower values of diffusion index 5 make the protons propagated spectra be harder for the same injection 
properties, resulting in the need for a softer proton injection index at high energies to reproduce the data, as shown 
in Table m where we show our resuhs for 5 = 0.33 ("KOL4-20"), S = 0.5 ("KRA4-20") and the intermediate value 
6 = 0.4 ("RUN4-20"). This in turn produces differences in the 7r° fluxes at the highest energies. 

Unlike protons, electron propagation at energies above 5 GeV is significantly affected by the energy loss time-scale 
and, since the ISRF and B-ficld model are kept fixed, the ICS and the higher part of the bremsstrahlung spectrum 
are not largely affected. In the very high energy part of the ICS spectrum we see a hardening for the models with 
greater S. That hardening is due to the fact that for greater S the higher energy diffuse faster out of the Galactic 
disk compared to lower energy e^, reaching the higher latitudes where we observe them through their ICS. In the 
lower energy part of the spectrum, bremsstrahlung varies significantly among the models since very different Alfven 
velocities and 77 values are used in those models in order to fit the CR data. While the overall fit of the 7-ray spectra 
is not affected much due to opposite effects of changing the value of S on the bremsstrahlung and the tt^ spectra 
below a few GeV, the relative ratio of bremsstrahlung to 7r° flux among the models with different S changes by up to 
a factor of two (at Ej ~ 0.5 GeV). Since both 7r° and bremsstrahlung are morphologically correlated to the gasses, 
discriminating among those components is very difficult. On the other hand, the predicted p spectrum favors larger 
values of S within our parameter search region. 

In Fig. Owe show the effect of varying the radial scale for the diffusion coefficient r^. Decreasing the value of 
results in lower values for the diffusion coefficient towards the Galactic center relative to the Sun's position, which 
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FIG. 3: Separate SNR dotted and pulsar dashed sources contributions to the diffuse 7-ray spectra for the Reference Astrophysical 
model "KRA4-20" . Sky windows are as in Fig. [J] 



forces the and p produced by sources closer to the Galactic center to spend a greater time close to the disk than 
those produced by sources close to the Sun. After refitting the diffusion coefficient normalization Dq (see Table [l| to 
the CR nuclear data, the net change in the fluxes is negligible. We also find that the difference between the radial 
independent case and the case with rd = 20 kpc is negligible. Yet the quality of the fit of the predicted p to the 
PAMELA p is affected by changing the scale rd and disfavors the smaller values of as is shown in Table HI 

In Fig. |6] we show the effect of varying the diffusion vertical scale Zd, that is correlated to the height of the diffusion 
zone |30|. Since smaller values of Zd yield a greater diffusion coefficient above the galactic plane, an overall rcscaling 
of the diffusion normalization Dq is necessary to fix the secondary to primary ratio (see Table |T|. Concerning protons, 
for which energy losses are less significant, changing Zd does not affect much their spectrum while also keeps the tt" 
spectrum and flux unchanged. Bremsstrahlung emission is also weakly affected by the changes in Zd-, because it is 
correlated morphologically to the gas distribution which is concentrated close to the galactic disk. Finally the ICS 
spectrum is mainly affected by the actual distribution of electrons being confined within thinner (thicker) diffusion 
zones resulting in lower (higher) total IC flux. For the case = 1 this results in a poor fit to the gamma-ray spectra. 
Thus thin diffusion zone models such as those that have been suggested by |102| , in order to give low p fluxes from 
Kaluza-Klein DM annihilation models, while simultaneously explaining the Icptonic excesses observed by ATIC and 
PAMELA, are in tension with the combination of CR and 7-ray spectra. 

In Fig. [7] we consider the effects of convective winds in the Galaxy. Convection introduces a new time scale into 
the propagation of CRs which mainly affects the protons, since energy losses still dominate electron propagation. 
After we refit to the B/C and proton fluxes, the diffusion properties are strongly affected (see Tab. |l| resulting in a 
significantly altered ICS and bremsstrahlung components. This is most evident at the low energy part of the spectra, 
where convection is more important relative to ICS and synchrotron losses. We find that high convection models are 
not favored by 7-ray data in the middle latitude region. In fact, we also find that low energy positron and electron 
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FIG. 4: Gamma-ray spectra predicted in models with dilTerent S. Plots refer to the different sky regions of our study, dotted 
lines: S — 0.5, dashed lines: S = 0.4, dashed-dotted lines: 5 = 0.33. For all cases = 4 kpc and rd = 20 kpc. 



fluxes arc in tension with PAMELA positron fraction data. 



C. Rigidity break in injection or diffusion 

Recently PAMELA jl03l | has observed a break at rigidity R ~ 230 GV, at both the proton and He spectra [STf . 
suggesting a hardening of the CR spectra at high energies. The harder spectral power-law at high rigidities is confirmed 
by the CREAM data as well [11]. Our combined fit of PAMELA and CREAM data leads to a break rigidity in our 
models at i?g 2 ~ 300 GV. 

One possible explanation for the observed rigidity break, is that the same break ori gina tes at the CR acceleration 
sites, at the SNRs shocks. Such a scen ario has b een suggested by studies of SNRs |104| | and from diffusive shock 
acceleration semi-analytical calculations jl05l4l09| . The pressure on accelerated particles arou nd th e shock leads to 
the formation of a precursor jl04l | where the upstream fluid is slowed down and compressed |l04| . For diffusively 
accelerated particles moving with respect to the shock, and thus between regions of different pressure, their gained 
energy depends on the " compression ratio" . On average, the higher energy particles which have larger diffusion 
lengths will probe the entire (or a greater part of the) precursor than the lower energy particl es, le ading to a concave 
shape spectrum. Thus the highest energy particles will "feel the total compression ratio" |l04| which (from first 
order F ermi accelerat i on) will result in the spectrum being harder than E~'^ at high energies and softer at low 
energies [lOl [TM [lOi, [ll3 . 

Another possible explanation is that at ~ 230 GeV we observe the emergence of a population of galactic sources 
(SNRs) that accelerate CRs with a resulting harder injection index. As long as this second SNR population is 
common enough in the Galaxy, and with a similar distribution as that injecting the softer CR spectra in the ISM, 
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FIG. 5: Gamma-ray fluxes for models in which we vary the diffusion radial scale rd. Plots refer to the different sky regions of 
our study, dotted lines: = 5 kpc, dashed lines: rd = 10 kpc, dashed-dotted lines: rd = 20 kpc. For all cases 5 = 0.5 and 
Zd = 4: kpc. 



both possibilities can be modeled in the same way with the DRAGON code, i.e. with the injection of CRs given by 
Eq. We will refer to this scenario as "scenario A", under which the 7-ray spectra of Fig. [^M] were produced. 

A third possibility is that instead we observe a change in the turbulence power spectrum of the ISM. The properties 
of the interstellar magnetic turbulence can be indirectly inferred from CR measurements. Before the PAMELA data, 
CR spectra were not measured accurately enough to exclude any break in the diffusion coefficient rigidity index. 
Coincidentally the needed change in the diffusion index A6 ~ 0.17 is the same as when considering the transition 
from a Kraichnan type turbulence at low i? to a Kolmogorov type at high R. We pursue this scenario ("scenario B") 
and after refitting the pulsar normalization flux to the total Fermi + e~ spectra, we calculate the 7-ray spectra, 
for I h \> 10°. In Fig. |S] we present the difference in the 7-ray flux between scenarios A and B, normalized to the 
fiux from scenario A, where we used the global diffusion model "KRA4-20" for our calculations. As it can be seen in 
Fig. m (left) the maximal difference between scenarios A and B, is up to 0(0. 1) in both ICS and 7r° spectra at energies 
^ 100 GeV, but of opposite sign, resulting in a difference in the total diffuse galactic component to be less than 
O(10~^) (Fig. [SI right) over the whole considered spectrum, including E^ ^ 100 GeV, in all three regions of interest. 
Such differences are too small to be probed by 7-rays at such high energies, because the extragalactic background 
fiux becomes more important, while also CR contamination modeling uncertainties, low statistics and the possible 
contribution from Dark Matte r ad d to the total uncertainty. 

Thus as also suggested by |lll| the best way to discriminate between scenarios A and B is through the p flux, 
where scenario A would give a soft break only at i? ~ 10 GV due to the break of the p and He spectra at ~ 230 GV, 
while scenario B would also give a harder break at ~ 230 GV (with the same spec tral index change as in the p and 
He fluxes) from the diffusion of the secondary p in the ISM (see discussion in |lll| ). 

Finally we note that [ll2j |. have discussed the impacts of a smooth hardening in the power-law of CRs and its 
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FIG. 6: Gamma-ray fluxes for models with different diffusion scale Zd- Plots refer to the different sky regions of our study. 
dotted lines: Zd = 1 kpc, dashed lines: Zd = i kpc, dashed-dotted Unes: Zd = 10 kpc. For all cases 5 — 0.5 and rd = 20 kpc. 



impact on high energy 7-rays and ps. After refitting the propagation parameters to the current wealth of data (which 
includes the most recent PAMELA data [131), and including the ICS and bremsstrahlung com pone nts, our differences 
in the total 7-ray spectra (shown in [5] left) turn out to be significantly smaller than those of jll2l |. 



D. Significance of the ISM Gasses 

As discussed in sectio n III Dl conventional models describing the HI and H2 interstellar gas distributions have been 
updated by the work of . l73l|. of which we show the profiles in Fig. [HI In Fig. 1101 using the assumptions of [70j for 
the HI, of [7^ for the H2 and of [t^ for the HII ISM gas, we show what are the separate contributions to the 7r° and 
bremsstrahlung 7-ray fluxes from each of those gasses. As it is clear, the HII contributes on average up to ~ 10% to 
the component fluxes at all energies and latitude regions that we show. Thus, changing the assumptions on the HI 
and H2 models can have a significant effect on both the propagation parameters shown in Tables HI and HII as has been 
sugg ested also in [H^l , as well as the 7-ray spectra. On the contrary the uncertainties on assumptions on HII (such as 
jll3| ) could at most result in a few % change in the 7-ray flux, well below the uncertainties from either one of the 
other two ISM gasses. Also note that since the protons suffer small energy losses their equilibrium spectrum is almost 
the same in the entire propagation region. Thus after changing the ISM assumptions and refitting the injection and 
propagation properties we also have similar CR protons density profiles in the Galaxy. Thus the changes that we 
observe in the 7-ray tt" fiuxes are very tightly correlated to the gas (target) distributions. In Fig. [11] we plot the 
steady state CR proton differential fiux profiles (density profiles) at i5 = 10 GeV (note that for higher energies the 
differences between the profiles, for given z or r, are even smaller). 

For the 7-ray spectra, the models of (tII ([zl]) for HI (H2) give -compared to the mean values of [t^, [t^- an 
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FIG. 7: Effects of convection on 7-ray spectra compared to the reference model. Plots refer to the different sky regions of our 
study, dotted lines: reference model, dashed lines: "CON4-20" model. 




FIG. 8: Absolute difference in fluxes between scenarios A and B (see text) normalized to the flux of scenario A. Left panel: tt" 
{dotted lines) and ICS {dashed- dotted lines). Red: 10° <| b \< 20°, Green: 20° <| b \< 60° and Blue: \ b \> 60°. Notice that 
the normalized differences for the tt" diffuse component have almost identical spectra. Right panel: Absolute difference in the 
total diffuse fluxes between scenarios A and B (including bremsstrahlung). Colors as in left-hand panel plot. The difference in 
the total diffuse are up to 2% (for Ej < 300 GeV), due to cancellation between vr" and ICS differences. 
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FIG. 9: Large scale density distributions of atomic (top) and molecular (bottom) hydrogen in the Galaxy vs r for z = (left); 
vs 2 for r = r© (right). For HI, "NS" refers to "GB" refers to [7l[ and "DL" to 0]- For H2, "Bronf refers to [tJ, "NS 
high" , "NS mean" and "NS low" refer to [t^ using, respectively, high, mean and low values of the midplane density. 



increased 7r° and bremsstrahlung flux by ~ 50% in the entire range of the spectra and at latitudes above 10° as shown 
in Fig. [T^l We note that 7-ray data favor the models [t^] and [t^ for the HI and H2 gas distributions, while the 
models of [tiI, [t^] and [zl] are marginally disfavored, from observations at the highest | 6 |> 60°, or intermediate 
10° <| 6 |< 20° latitudes (see Table |T|. This is mainly due to the higher density of the local H2 local, as shown in 
Fig. |9] (lower row). 

That can also be understood from Fi g. [1^ (lower right panel) where we plot the relative difference in the 7r° 7- 
ray fluxes between the H2 gas model of [TJI and the reference H2 model. The differences in the 7r° fluxes between 
these models are of 0(1). We also note that in the original parametrization of [zl], the assumed Xco factor = 
N(H2)/W(CO) was taken to be constant, equal to (2.8 ± 0.4) x 10^° cm-^K-^km'^s. This is a factor of 2 higher 
than the Xco fact or of Eq. (TUl that is in better agreement with recent 7-ray analysis of the ISM in the outer part 
of the Galaxy [ll4l |. Had we used the constant value for the Xco factor, our 7r° and bremsstrahlung fluxes would be 
enhanced, bringing them in more tension with the 7-ray data. Also we alternatively use the la higher and la lower 
values for the H2 density proflle in r of [t^, which we plot also in Fig. [T^ and CTbottom) . using however the mean 
values of [t^I for the HI gas. We note that both the la higher and the la lower cases are in tension with the data as 
well. 

From Fi g. [1 21 (lower right panel) it can be also seen that switching from the older parametrization [7l|, [T^I , to the 
newer one [701 ^'^^ the HI gas has much smaller effect compared to switching between H2 models. The reason for these 
much smaller differences in HI is that the steeper decrease with distance from the Galactic plane present in [t^I (see 
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FIG. 10: The observed 7r''(reci) and bremsstrahlunR (green) components of the diffuse 7-ray flux, from the three major ISM 
gas components. Dotted lines: from 112 following [73| (mean values), dashed lines: from HI following [70[, dashed-dotted lines: 
from HII following [zl, solid lines: total flux. Upper left: 10° <| b \< 20°, upper right: 20° <| h \< 60°, lower left: \ b \> 60° 
averaging over all latitudes. Lower right: flux ratios between the components for | 6 |> 60° to 10° <| 6 |< 20°. 



Fig [5] top right panel) as o pposed to [t^I , compensates for its suggested higher density on the disk [7l| (see Fig. [3] top 
left panel) . Recently |ll5l | has also shown the importance of the uncertainties of the gas models and especially that 
of the Xco f actor , suggesting as we do that the uncertainties in the H2 gas distribution are the greatest (see relevant 
discussion of |115| ). 

The increased number of target nuclei in the gas models [tiI, [tI, [tI] , results in the need of a faster escape of CR 
nuclei from the Galaxy, in order not to overproduce secondaries. This will then increase the diffusion coefficient 
normalization to keep the same B/C flux ratio. Therefore would propagate to larger distances from the Galactic 
disk, resulting in a 10% increase of the observed IC flux, as well as to an effect on the tt" and bremsstrahlung 
components. 

A possible way to reconcile an increased number of target nuclei in the ISM gas would be to decrease the thickness 
of the diffusion halo. That is shown in Tabic HII where we show the effects of interplaying the gas distribution and the 
thickness of the diffusion halo. Our model "KRAl-20 Bronf ' represents our thinnest diffusion halo with the highest 
ISM gas assumption, which can not be considered in tension with any data. Yet the reverse case, of a thick halo with 
a low ISM gas ("KRAlO-20 NS low") tends to undcr-prcdict the 7-ray spectra at low latitudes. Finally, the extreme 
cases of a very thin diffusion halo and a low ISM gas assumption ("KRAl-20 NS low"), or a high ISM gas in a thick 
diffusion halo ("KRAlO-20 Bronf), while can still predict the p flux in good agreement to the data, are in systematic 
tension with the 7-ray fluxes by either systematically under-predicting them or systematically over-predicting them. 
Thus a combined analysis of CRs and 7-rays as is ours can probe the uncertainties in the large scale gas distributions, 
which when using only CR information, are large, especially for the H2 gas. 
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FIG. 11: CR proton flux profile at _E = 10 GeV (at steady state) as a function of r for given z (left) and as a function of z for 
given r (right). In "Bronf we use the model of [tJ for H2 and the model of [tH. It^] for HI gas distributions. In "NS high" , 
"NS mean" and "NS low" we use for H2 gas distribution, respectively, high, mean and low values of midplane density of [73 ]. 
however they share the same HI gas distribution of [70l |. 



E. SNR Distribution 



Since in our study we use the 7-ray data above | ^ |> 10° we are not very sensitive to the SNR source distribution. 
In particular, we are poorly sensitive to the inner few kpc, that are also hardly probed by direct observations of single 
sources. To study the significance of the source distribution in the inner part, we use three different source profiles 
as shown in Fig.[T31 The Ferriere et al. [s^ profile ("Source A"), our reference assumption, is recovered from Type 1 
and Type 11 supernovae distribution models and is defined by 

fs{r,z) = o.l38e-('-'^«)/4-5-l^l/°-325 

+ (0.79e-(^/O-2i2)^ o.21e-(^/°-636)^) x OMSe-^^'-^-e)/^-^' r > 3.7 kpc 

+ (0.79e-(^/°-2i2)^ + 0.21e-(^/0-63e)^^ ^ 3.349e-('-3-^)'/2.i= ^ < 3 7 ^^4) 

The relative normalization of the t wo p opulations is based on the averaged occurrence frequencies of SNe Type 1 
and Type II in other galaxies [llTl Ill8l | , while the spatial profiles are based on the assumption that Type 1 have a 
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FIG. 12: Effects of changing the models that describe HI and H2 gas distributions on the 7-ray spectra in the sky regions of our 
study, dotted lines: HI from [tH. It^ and H2 from [tJI- solid lines: HI from [70| and H2 from [T^l, dashed lines: HI from [tO] 
and H2 from [tJ increased by 1 a, dashed-dotted lines: HI from [t^] and H2 from [t^], diminished by 1 a. In the lower right 
panel we show the relative difference between the tt" components predicted by various gas assumptions and the one predicted 
by our reference model. Red: 10° <| 6 |< 20°, green: 20° <| 6 |< 60° and blue: | 6 |> 60°. The dotted lines are practically 
overlapping due to the H2 models of [T^I and [TJI having a similar vertical scaling. 
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Red: parametrization of [3y| given in Eg. 1141 green: discribed by Eg. 1151 [116t | and blue: given by 
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distribution similar to that of old disk stars |119{ , while Type II (which are the most frequent) are tightly correlated 
to the arms. Especially in the inner 3.7 kpc the profile of the Type II is correlated to our Galaxy's pulsar distribution 
and thus sensitive to selection effects. 

Also the distribution given in Eq.[T5]lias been used extensively in GALPROP [60llll6j as a conventional distribution 
for SNRs. 

This distribution parameterizations ( "Source B" ) comes from actual observations of galac tic SNe |l20l | with its detailed 
values being selected to agree with the diffuse fluxes of analyzed EGRET 7-ray data jll6l Il2lj . As a result, this 
parametrization has to be taken with some care, in particular with respect to its prediction on the inner parts of the 
Galaxy. 

Since both parametrizations of [36| and jll6l | are least predictive towards the inner regions of the Galaxy, we study 
also a third parametrization described by "Source C" : 

/p(r, z) = (0.078 + 2.57e-(''/'"0)') e-l^l/°-2. (16) 

As is also shown in Fig. I13[ this parametrization gives an almost constant radial distribution in the inner 3 kpc, 
while its averaged distribution at r > 5 kpc is similar to that of "Source A" and "Source B" , with its vertical behavior 
being the same as that of "Source B" . Since it is expected that Type II SNe are correlated to the spiral arms, and 
since the inner few kpc are populated by old stars, we use "Source C" as a probe of the maximal effect that the 
uncertainties in the SNRs distribution could have on the diffuse 7-ray data analysis, rather than as an optimal model 
for SNRs. Yet as it can be seen from Table U and Fig. [TH where we plot the total 7-ray fluxes for the three SNR 
models, "Source C" indeed provides a better fit to the 7-ray data at low latitudes, that probe the inner parts of the 
Galaxy. That could be an indication that either "Source A" and "Source B" , being connected to observations, indeed 
under-predict the "recent" SNR density towards the inner part of the Galaxy, or simply, as suggested in section HVAl 
that unresolved point sources are present. Also changing the gas assumption only in the inner few kpc could have an 
effect. To probe such uncertainties a study including lower latitude regions toward the GC would be very well suited. 



F. pp-Collision 7-ray spectra 

In Fig.[Tn]we show the 7r° component of the diffuse 7-ray fiux at 10° < | | < 20° using three different parametrizations 
for the 7-ray spectra produced by pp collisions. We clarify that we call these spectra "tt" spectra" since their 
contribution to the 7-ray spectra is the dominant one, in all parameterizations. However also the 7-rays from the 
decay channels of other produced mesons such as K^, K°, 77, D\D° are taken into account. 

As a reference parametrization we use that of Kamae et al. [l22l Il24{ that was derived for the cross-sections of 
diffractive, non-diffractive and excitation of resonance processes, base d on simulation and experimental data o n pp 



collisions. We compare that parametrization to that from Kelner et al. [123l | that was based on running SIB YLL [125 1 
simulations of pp collisions and also to the 7-ray spectrum from our own P ythia (version 6.4) simulations jl26l |. For 
the Ka mae et al. parametrization we u se the updated tables information jl24| relev ant to those of tables 2 and 3 
of |l22| that ar e us ed in eq. 5-14 of jl22l |. For the Kelner et al. parametrization jl23l | we used the information given 
in eq. 58-61 of |l23l |. while in our Pythia simulations we run pp collisions with center of mass energy from 2.33 GeV 
up to 7 TeV, with subsequent decay of all mesons and including final state radiation. We keep the information for 
the 3D momenta of the final stable particles, which we re-boost to the proper observer frame (where a CR p hits 
a practically stable ISM p). We see that the 7-ray spectra, normaliz ed a t 10 GeV, agree well from 100 GeV down 
to energies of 1 GeV where it is expected that the simulations from |l23l | and Pythia would be no longer reliable. 
Thus it is safe to say that uncertain ties in the 7-ray spectra produced by pp collisions that could be due to missing 
processes in the parametrization of jl22j |. are too small to have a strong impact on the constraints imposed on the 
ISM properties that we have described using the combined CR and 7-ray spectra. 



V. CONCLUSIONS AND DISCUSSION 



The results of the analysis show that combining the local CR measurements, which are rather powerful to constrain 
the local averaged properties of propagation and source/gas distribution, with the diffuse gamma-ray fluxes at inter- 
mediate and high latitudes in the energy range currently covered by Fermi, can be useful in giving some constraints 
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FIG. 14: Gamma-ray fluxes assuming different models for the SNR distribution. Plots refer to the different sky regions of our 
study. For each SNR model we r efit the diffusion parameters, as is shown in Table ID Dotted lines: Source A of Eq. [14] [s^ . 
dashed lines: Source B of Eg. ITS] [116|] . dashed-dotted lines: Source C of Eg. \TS\ 




FIG. 15: The diffuse tt" component of the flux using the Kamae et al. [l22| para metri zation for pp collisions vs that from [l23l ] 
and Pythia runs. Red solid: Kamae et al. parametrization. Blue dotted: using [123| and blue dashed: it s diff erence from the 
[122| [. Green dashed-dotted: using Pythia, and green dashed- dotted- dotted- dotted its difference from the [122| 1. We normalize 
the tt" fluxes at 10 GeV. 
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also on the global galactic parameters. We have tested a large subset of the benchmark models selected to depict 
rather diverse settings for modeling propagation and the ISM. While many of these models provide a good fit to both 
the CR and the 7-ray data, still there are a few scenarios that our analysis can disfavor. 

In particular, we have studied rather extreme limits on CR diffusion galactic profiles, ranging from essentially 
constant diffusion coefficient everywhere in the Galaxy, down to exponentially suppressed vertical and radial profiles 
with I kpc scale hight and 5 kpc radial scale. The combined fit of CRs and 7-rays suggests a slight preference for 
thicker diffusion zones, while there is a weak dependence on the variation of the diffusion coefficient in the radial 
direction, which is however better probed and constrained by the antiproton spectrum. While this result is not 
conclusive, it suggests a trend that better statistics and smaller systematics on the 7-ray spectra (soon to come) will 
lead in testing specific models for the position dependence of the diffusion coefficie nt D (r, z) of CRs in the Galaxy. 
The high statistics measurements of the local flux of radioactive isotopes by AMS-02 |l27l | , placed on the International 
Space Station, will add further information on the vertical thickness of the diffusion region, possibly allowing to break 
the degeneracy between thicker regions of emissivity populated by CRs diffusing out of the galactic disk and exotic 
sources with an intrinsically thicker scale height, such as from DM. 

Moreover, the current diffuse 7-ray spectra can discriminate (and even constrain) profiles for the ISM gasses. Before 
Fermi-hAT 7-ray data, in order to place constraints on the ISM properties, the secondary to primary CR spectra were 
used, with the best measured data sets coming from B/C and p/p. Yet, with only the CR data as a handle, a variation 
of the large scale gasses distributions could almost always be compensated by changing the diffusion properties (mainly 
the normalization of the diffusion coefficient). By exploiting the 7-ray diffuse fluxes above | & |> 10° and combining 
them with the CR data, we have shown that we can actually break the degeneracy between diffusion and ISM gas 
distribution. In fact, thanks to the expected improvement both in statistics and systematics errors of the 7-ray data 
from Fermi-LAT, and even more with the CR spectral data up to Fe from 0.1 GeV/n to at least 100 GeV/n from 
AMS-02, we can be optimistic in further constraining the properties of the ISM gas distributions, within the next few 
years. 

On the other hand, we also find that CR and 7-ray data do not constrain strongly the diffusion spectral index 5 
within the range we considered. 

Furthermore, we have discussed the implications from the recently found rigidity break in the protons and He CR 
spectra [131 (confirmed also by H^). We have addressed the possibility of discriminating whether the break is in the 
injection spectrum (connected to either acceleration effects in the sources, or to the presence of an extra population of 
primary sources injecting CRs with harder spectra) or in the energy dependence of the diffusion coefficient. We have 
found that the galactic diffuse 7-rays cannot be used to this aim, neither with the current nor with the near futu re 
projected accuracy of the spectra, leaving this task to other observables, such as antiprotons as suggested by [lll| . 

As a final remark, we have shown that our analysis is robust with respect to uncertainties in the parameterization 
of the 7-ray spectra produced in pp collisions. 

Having achieved with these results a better understanding of the contributions of the astrophysical components to 
the diffuse 7-rays, a natural application will be to place limits on a possible exotic contribution to the high latitude 
7-ray flux. Our forthcoming analysis devoted to the search for a DM signal will be of particular relevance in this 
respect. 
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Appendix A: Synchrotron losses of CR electrons 

Cosmic Ray electrons as they propagate in the Galaxy, lose energy via synchrotron and bremsstrahlung radiation, 
and by up-scattering low energy photons to higher energies (as X-rays or 7-rays) or by ionization losses inside gasses. 
Adiabatic losses can also be important inside expanding SNRs but do not matter when CRs are considered free to 



* Smaller scale features are much better probed by synchrotron data, as for example has recently been done in |128| . 



24 




FIG. 16: Ratio of the energy loss rate of 1 GeV{top left), 10 GeV {top right) and 100 GeV (bottom) electrons due to synchrotron 
radiation to the total energy loss rate. 

diffusively propagate in the interstellar medium (and not inside some expanding volume of matter). Above energies 
of several GeVs, energy losses via inverse Compton scattering and synchrotron radiation become most important, 
while at lower energies and denser ISM environments bremsstrahlung and ionization can dominate. 

As we have described in section IIIBl we expect a correlation between the large scale profile of the galactic B-fields 
and that of CR diffusion. Yet, in this work, for simplicity we have chosen to keep a specific parametrization of 
the B-field as given in eq. [5] (following WMAP [EH), while changing the diffusion coefficient's profile described by 
eq. [9] Such a simplification is possible since the B-field assumptions in our code can influence only the synchrotron 
emissivity calculations, which are not presented here and are left for future work, and the synchrotron energy losses 
of CR electrons. The latter are subdominant at all energies and positions of the Galaxy with respect to the energy 
losses due to up-scattering CMB, infrared and optical photons, or bremsstrahlung radiation, by at least a factor of 3 
for typical values of galactic magnetic fields. For the magnetic field that we used in this paper, we show in Fig. [16] 
the ratio of synchrotron energy losses rate to the total energy losses rate for electrons oi E = 1, 10 and 100 GeV. 

As can be seen, the synchrotron energy losses are in no part of the Galaxy, dominant, accounting at most up to 
22% of the total electron energy losses at 100 GeV, locally. The synchrotron radiation losses are more important at 
very high energies, where the ICS losses due to optical and infrared photons up-scattering become less efficient as the 
Klein-Nishina cross-section decreases away from the Thomson cross-section value. Also, far away from the galactic 
disk, synchrotron radiation losses (which scale with the square of the magnetic field strength) drop because of the 
exponential decrease of the B-field (eq. [8]), while the energy losses due to CMB up-scattering remain the same. 

Therefore, we can treat the effects of the magnetic fields on diffusion and on energy losses to a good approximation 
separately. 

Appendix B: Impact of a 2D vs a 3D ISM gas distribution on the diffuse 7-ray spectra 

The Fermi-h AT ^ 1° angular resolution (for energies above few GeV), allows to trace in some detail the morphology 
of the 7-ray emissivity associated to the gas in the Galaxy. Thus one can use a 3-dimensional ISM gas distribution 
to compare to the 7-ray data. In this paper we have instead used 2-dimensional spatially smoothed gas distributions, 
which does not account for any of the small-scale features in the 7-ray maps. 
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FIG. 17: Gamma-ray spectra for the 3 regions of interest using the 3D formahsm in DRAGON. 



Yet, our goal has not been to study or interpret these structures, but rather the larger scale properties of the Galaxy, 
that are incorporated in the 7-ray spectra measured in the very wide angular windows that we use. To illustrate the 
minimal impact on our analysis of having averaged out the small-scale features, we consider our reference propagation 
model and calculate in the three regions: 0° < I < 360°, 10° <| 6 |< 20° / 20° <| 6 |< 60° / I ^ l> 60°, the diffuse 
7-ray spectra using a 3D model for the HI and H2 gas, equivalent to the 2D reference model; results are shown in 
Fig. 1171 and should be compared to those in Fig. [2] obtained under the same configuration for cosmic-ray propagation 
but with the 2D gas model. 

The differences in the tt" and the bremsstrahlung (relevant for ISM gas) arc at the 10% level, in the three parts of 
the sky that we study. This is most clearly seen in Fig. [THl where we show the normalized differences in using the 2D 
and the 3D gas distributions. 

These results are well within the accuracy needed for our analysis, in which wc have studied models of HI and H2 
gas distributions that can differ by up to 30% in their prediction of the 7r° and bremsstrahlung diffuse fluxes (Fig. [10] 
bottom right). 

The total galactic diffuse model predictions between the 2D and the 3D cases are actually less than 5%, as there 
arc compensations between the different components. Furthermore, given that in deriving our physical conclusions we 
use the analysis carried between the total 7-ray Fermi fluxes and the total 7-ray predicted fluxes, which include 
the EGBR and point sources, the impact of using a 2D ISM gas model is minimal. Finally using our 3D ISM gas 
model, we have checked that the contribution from the " Fermi Bubbles" [l^] (that we do not include explicitly) in the 
windows of 0° < Z < 360° is of roughly the same magnitude as the numerical uncertainties in calculating the 7-ray 
fluxes. 
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